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Abstract 



>. 

CN , We introduce a particle-number reprojection method in the shell model 

Monte Carlo that enables the calculation of observables for a series of nuclei 

q^ , using a Monte Carlo sampling for a single nucleus. The method is used 

ON 

to calculate nuclear level densities in the complete (pf + gg/ 2 )-shell using 

a good-sign Hamiltonian. Level densities of odd-^4 and odd-odd nuclei are 
Q ■ 6 6 

reliably extracted despite an additional sign problem. Both the mass and the 



T z dependence of the experimental level densities are well described without 



any adjustable parameters. The single-particle level density parameter is 
found to vary smoothly with mass. The odd-even staggering observed in the 
calculated backshift parameter follows the experimental data more closely 
than do empirical formulae. 
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The interacting shell model has successfully described a variety of nuclear properties. 
However, the size of the model space increases rapidly with the number of valence nucleons 
and/or orbits, and exact diagonalization of the nuclear Hamiltonian in a full major shell 
is limited to nuclei with A < 50 fl]^. The development of quantum Monte Carlo meth- 
ods for the nuclear shell model allowed realistic calculations of finite- and zero-temperature 
observables in model spaces that are much larger than those treated by conventional diago- 
nalization techniques @,[|. 

The Monte Carlo method was successfully adapted to the microscopic calculations of 
nuclear level densities 0. Accurate level densities are needed for estimating nuclear reaction 
rates, e.g., neutron and proton capture rates. The nucleosynthesis of many of the heavy 
elements proceeds by radiative capture of neutrons (s and r processes) or protons (rp process) 
in competition with beta decay 0,0]. Theoretical approaches to level densities are often 
based on the Fermi gas model, i.e., the Bethe formula ||, which describes the many-particle 
level density in terms of the single-particle level density parameter a. Shell corrections 
and two-body correlations are taken into account empirically by defining a fictitious ground 
state energy. In the backshifted Bethe formula (BBF) the ground state energy is shifted by 
an amount A. This formula describes well the experimental level densities of many nuclei 
if both a and A are fitted for each individual nucleus ||. While these parameters have 
been discussed in terms of their global systematics, it is difficult to predict their values for 
particular nuclei. 

The nuclear shell model offers an attractive framework to calculate level densities, but 
the model space required to calculate level densities at excitation energies in the neutron 
resonance regime is usually too large for conventional diagonalization methods. We have 
recently developed a method to calculate exact level densities using the shell model 
Monte Carlo (SMMC) approach, and applied it to calculate the level densities of even- 
even nuclei from iron to germanium ||10| . Fermionic Monte Carlo methods are usually 
hampered by the so-called sign problem, which causes a breakdown of the method at low 
temperatures. A practical solution in the nuclear case was developed in Ref. @] but the 



associated extrapolation errors are too large for reliable calculations of level densities. In 
Ref. |5[] the sign problem was overcome by constructing a good-sign interaction in the (pf + 
<79/2)-shell that includes correctly the dominating collective components of realistic effective 



interactions |iT|. The SMMC level densities are well-fitted by the BBF, and both a and A 
can be extracted from the microscopic calculations. 

The SMMC approach is computationally intensive. In particular, the SMMC level den- 
sities require calculations of the thermal energy at all temperatures. The weight function 
used in the random walk is temperature dependent, and a new Monte Carlo sampling is 
required at each temperature. Since this procedure has to be repeated for each nucleus, the 
calculations are time-consuming. In this paper we describe a particle-reprojection method 
that allows us to calculate observables for a series of nuclei using Monte Carlo sampling 
for one nucleus only. The random walk is done with a weight function proportional to the 
partition function of a given even-even nucleus (which is positive definite for a good-sign in- 
teraction), and the thermal observables are then calculated for several nuclei by reprojection 
on different particle numbers (both even and odd). This method allows significantly more 
economical calculations of level densities. We apply the method in the full (pf + p 9 / 2 )-shell 
to study the systematics of a and A for even-even, odd- A and odd-odd manganese, iron and 
cobalt nuclei. A direct comparison with both experimental data and empirical formulae is 
presented. The agreement with the data is remarkably good with no adjustable parameters 
in the microscopic calculations. Furthermore, we find that the SMMC values follow the data 
more closely than do the empirical values. 

The Monte Carlo method is based on the Hubbard-Stratonovich representation of the 
imaginary-time propagator, e~^ H = J D[a]G(a)U (T , where G(a) is a Gaussian weight and U a 
is the propagator of non-interacting nucleons moving in fluctuating auxiliary fields a which 
depend on imaginary time. The canonical thermal expectation value of an observable O is 
given by (O)^ = J D[(r]G(a)Tr^(OU (r )/ J D[a]G(a)Tr^U (T , where Tr_4 denotes a trace in the 
subspace of a fixed number of particles A. In actual calculations we project on both neutron 
number N and proton number Z, and in the following A will denote (N, Z). We rewrite 



(0) A = (Tr A (OU a )/Tr A U a ) w , (1) 

where we have introduced the notation 

_ JD[a}W(a)X a 
{X ° )w = JD[a]W(a) ' (2) 

and W(a) = G{a) r Ti A U a . For an even number of particles with a good-sign interaction, 

W(a) is positive definite. In the Monte Carlo method we choose M samples (each denoted 

by (Tfc) according to the weight function W(a), and estimate (X a )w ~ J2kX ak /M. 

We assume that the Monte Carlo sampling is done for a nucleus with particle number A, 

and consider the ratio Z A //Z A between the partition function of a nucleus with A' particles 

and the partition function of the original nucleus. In the notation of Eq. (Q) 

ZAP) _Tr A ,e-P H /Tr A ,U a > 



Z A ((3) Tiue-ZMT \Tv A Uj w - (3) 



Similarly, the expectation value of an observable O for the nucleus with A' particles can be 
calculated from 

\ Q IA' = / Tr A ,[/ CT \ • ( 4 ) 

\Tr^C/ CT / w 

The Monte Carlo walk is carried out by projection on a fixed A, and Eqs. (^) and (^) are 
then used to calculate the partition functions and observables for a family of nuclei with 
A'^A. 

We applied the method to nuclei in the (pf + g 9 / 2 )-shell, using the Hamiltonian of 
Ref. ||]. The single-particle energies are computed in a central Woods-Saxon poten- 
tial V(r) plus spin-orbit interaction, while the two-body interaction includes a monopole 
isovector pairing of strength go plus a separable surface-peaked interaction Jl2| v(r,r') = 
—x{dV/dr)(dV/dr')5(r — r). The surface-peaked interaction is expanded into multipoles 
and only the quadrupole, octupole and hexadecupole terms are kept. The strength \ is deter- 
mined self-consistently and renormalized. The strength of the pairing interaction go w 0.2 is 
determined from experimental odd-even mass differences. Both the pairing and the surface- 
peaked interactions are attractive and lead to a good-sign Hamiltonian. A repulsive isospin- 



dependent interaction leads to a sign problem, and was included perturbatively in recent 
level density calculations in sci-shell |13] and p/-shell |TJ]] nuclei. 



In the particle-reprojection method described above we have assumed that the Hamilto- 
nian H is independent of A. Suitable corrections should be made if some of the Hamiltonian 
parameters vary with A. Since \ depends only weakly on the mass number A (oc A^ 1 ^ 3 ), 
and the pairing strength go is constant through the shell, the largest variation is that of the 
single-particle energies of the orbit a, t a {A). To correct this variation we approximate the 
thermal energy of A' = (N', Z') particles by 

Ea'(P) « £kGA') - e a (A)](n a ) A , + (H) A , (5) 

a 

where H is the Hamiltonian for a nucleus with A particles. In calculating the energy for 
A' particles from (||), we used in the propagator {e~^ H ) the Hamiltonian H for nucleus A 
rather than A'. To minimize the error we reproject on nuclei with N' — Z' values close to 
N — Z (the Woods-Saxon potential depends on N — Z). In the applications below we have 
checked that the resulting error in the level density is negligible. 

We used the reprojection method to calculate the thermal energies versus (3 for 50_56 Mn, 
52 ~ 58 Fe, and 54_60 Co including odd- A and odd-odd nuclei. We sampled according to the 
even-even nucleus 56 Fe and reprojected onto 53_56 Mn, 54_58 Fe, and 54_60 Co, while the nuclei 
50_52 Mn and 52 > 53 Fe are reprojected from Monte Carlo sampling of the odd-odd N = Z 
nucleus 54 Co. The calculations are done for values of f3 between (3 = and 1 MeV -1 in 
steps of A/5 = 1/16, and between 1 and 2.5 in steps of A/? = 1/8. At each j3 we used about 
4000 independent samples. Reprojected energy calculations typically have larger statistical 
errors at larger values of (3. Therefore we also performed direct Monte Carlo runs (without 
reprojection) for the above nuclei at several values of (3 between 1.75 and 3.0 MeV -1 . For 
odd-A and odd-odd nuclei, a typical statistical error for the energy at /3 ~ 2.5 is ~ 0.5 
MeV, while for (3 ^ 3 the error is too large for the data to be useful. This is just another 
manifestation of the sign problem for nuclei with an odd number of protons and/or neutrons. 
Fortunately, because of the high degeneracy in the vicinity of the ground state of these nuclei, 



the thermal energy is already close to its asymptotic value. 

Fig. [I] shows the calculated SMMC thermal energies versus (3 for a series of cobalt 
isotopes. The effect of pairing on the thermal energies at low temperatures (i.e. large j3) is 
clearly seen in their uneven spacings. The inset of Fig. [I] shows the SMMC thermal energies 
(triangles with error bars) for 60 Co for the large values of /3 only. 

In calculating the level density versus excitation energy, it is important to get accurate 



values of the ground state energy. In Ref. |L0| we used, for even-even nuclei, a two-state 
model (0 + and 2 + ) to obtain a two-parameter fit to the thermal energy and (J 2 )- For odd- 
A and odd-odd nuclei this method is not useful since in general we do not know the spin 
of the ground state and first excited state. Moreover, these nuclei do not have a gap and 
often more than two levels contribute to the thermal energy at the lowest temperatures for 
which Monte Carlo calculations are still possible. We estimate the ground state energy of 
these nuclei by taking an average of the large-/? SMMC values for the thermal energy. The 
diamonds in the inset of Fig. [I] are such average values for 60 Co. We estimate the ground 
state energy of the odd- A and odd-odd nuclei to be reliable to about ~ 0.3 MeV. 

To calculate the level density we first find the partition function Zj\> by integrating the 
relation — din Z A >/df3 = E^. The level density is then given by 

p A , = (2n[3- 2 C A ,)- 1/2 e s * , (6) 

in terms of the canonical entropy S A > = /3E A i + In Z A i and the heat capacity C A i = 
-(3 2 dE A ,/df3. 

The level densities for the cobalt isotopes of Fig. [I] are shown in Fig. ^| as a function of 
excitation energy. These total level densities are fitted to 

p(E x ) « g^a-*(E x - A + t )-i e V«^-A> , ( 7 ) 

where t is a thermodynamic temperature defined by E x — A = at 2 — t and g = 2. Eq. 
(0) is a modified version of the BBF derived by Lang and Le Couteur |JL5|| . It differs from 
the usual BBF by the additional "temperature" term t in the pre-exponential factor, and 



provides a better fit to the calculated level density at lower excitation energies. The solid 
lines in Fig. ^| are the BBF level densities (|7|) fitted to the SMMC level densities in the range 
E x < 20 MeV. In general we obtain a good fit down to energies of ~ 1 MeV or smaller. The 
inset shows the low energy fit for 55 Co. The dashed line is a fit to the BBF without t; this 
approximation starts to diverge around 2 MeV due to the singularity of the pre-exponential 
factor {E x — A)~ 5//4 . Notice that the level density for an odd-odd cobalt (e.g., 54 Co) is higher 
than the level density of the subsequent odd-even cobalt (e.g., 55 Co) even though the latter 
has a larger mass. This is due to reduced pairing correlations in the odd-odd nucleus that 
lead to a smaller backshift A. 

We extracted the level density parameters a and A for the above nuclei by fitting Eq. 
(0) to the SMMC level densities. The results for a and A versus mass number A are shown 
in Fig. ^. The Monte Carlo results (solid squares) are compared with the experimental 
data (X's) quoted in Refs. || and flllj. The solid lines describe the results of the empirical 



formulae of Refs. |I7|] . The calculated values of a depend smoothly on the mass, unlike 
some of the empirical results, and in the case of the cobalt isotopes follow the data more 
closely. The staggering seen in the behavior of A versus A is a result of pairing effects. In 
the empirical formulae, A ~ for odd-even nuclei, is positive for even-even nuclei and is 
negative for odd-odd nuclei. We see that the present values of A follow rather closely the 
experimental values, and are in general in better agreement than the empirical values. The 
lower values of a (relative to the experimental values) for the odd-odd manganese nuclei 
are compensated by corresponding lower values of A, and thus do not cause significant 
discrepancies in the level densities for E x < 10 MeV. 

To demonstrate how the Monte Carlo results improve over the empirical formulae we 
show in Fig. |] the calculated level densities of three A = 55 nuclei (Mn, Fe and Co). 
According to the empirical formula, A ~ for odd-A nuclei, and the values of a are similar 
(since A is the same). The empirical formulae therefore predict similar level densities for 
these nuclei. However the SMMC level densities of these three nuclei are seen to be quite 
different from each other. Indeed we find that A is positive for 55 Co, close to zero for 



55 Fe and negative for 55 Mn. The experimental level densities (dashed lines) are in good 
agreement with the Monte Carlo calculations, suggesting a T z = (N — Z)/2 dependence of 
the level density, which is usually ignored in empirical formulae but is clearly observed in 
our microsocpic calculations. 

In conclusion, we have described a particle-reprojection method in the SMMC that en- 
ables the calculation of thermal properties for a series of nuclei using Monte Carlo sampling 
for a single nucleus. We applied the method to the calculation of level densities. 
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FIG. 1. The SMMC thermal energies versus (3 for 54_60 Co (symbols). Shown on the right are 
the extrapolated values of the ground-state energy. Inset: the SMMC thermal energies (triangles) 
at large f3 values for 60 Co. The diamonds are the energies obtained by averaging the large- f3 results 
above the corresponding (3. 
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FIG. 2. The SMMC level densities of the °Co isotopes. The solid lines describe a fit to 
the BBF (|7|). The top inset shows the level density of 55 Co at low excitation energies. The circles 
are the SMMC results, the solid line is a fit to (0), and the dashed line is a fit to (0) but without 
the t term. 
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FIG. 3. The single-particle level density parameter a (left column) and the backshift parameter 
A (right column) for Mn, Fe and Co isotopes. The solid squares (connected by dashed lines) are the 
results of fitting the calculated SMMC level densities to Eq. (0), while the X's are the experimental 
results. All the experimental values are from the compilations of (assuming rigid body moment 
of inertia), except for 58 Co and 59 Co where the values quoted in [l(| are used. The solid lines are 



the empirical formulae of [17]. 
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FIG. 4. The SMMC level densities of three A = 55 nuclei: 55 Mn (circles), 55 Fe (squares) and 
55 Co (diamonds). The dashed lines are the experimental level densities. 
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